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Abstract — Relay-based cooperative communication has become 
a research focus In recent years because it can achieve diversity 
gain In wireless networks. In existing works, network coding 
and two-path relay are adopted to deal with the Increase of 
network size and the half-duplex nature of relay, respectively. 
To further Improve bandwidth efficiency, we propose a novel 
cooperative transmission scheme which combines network coding 
and two-path relay together In a multi-source system. Due to 
the utilization of two-path relay, our proposed scheme achieves 
full-rate transmission. Adopting complex field network coding 
(CFNC) at both sources and relays ensures that symbols from 
different sources are allowed to be broadcast in the same time 
slot. We also adopt physical-layer network coding (PNC) at relay 
nodes to deal with the inter-relay Interference caused by the two- 
path relay. With careful process design, the ideal throughput 
of our scheme achieves by 1 symbol per source per time slot 
(sym/S/TS). Furthermore, the theoretical analysis provides a 
method to estimate the symbol error probability (SEP) and 
throughput in additive complex white Gaussian noise (AWGN) 
and Rayleigh fading channels. The simulation results verify the 
improvement achieved by the proposed scheme. 

Index Terms — complex field network coding; two-path relay; 
inter-relay interference; physical-layer network coding; through- 
put. 



I. Introduction 

Based on the fact that it is difficult to pack multiple anten- 
nas per terminal due to the limit of physical space, relay- 
based cooperative communication is especially eye-catching 
since it can extend coverage and achieve spatial diversity 
gain in wireless communication networks. However, due to 
the increase of network size and the half-duplex nature of 
relay, traditional relay schemes are bandwidth inefficient. 
Two methods were proposed to break through this bandwidth 
bottleneck: 

1) Network coding. Physical-layer network coding (PNC) 
was proposed in [ 1 ] to map superposition of electromagnetic 
signals to simple Galois field GF(2") additions of digital 
bit streams. In [2], Galois field network coding (GFNC) 
was applied to wireless communication networks to enhance 
bandwidth efficiency. Afterwards, Wang and Giannakis put 
forward a novel cooperative approach based on complex field 
network coding (CFNC) [3], which makes it possible for 
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Time slot n 
Fig. 1. Two-path relay scheme model 



sources to broadcast symbols in the same time slot (TS) on 
the same resource blocks (RBs) after symbol-level operations 
at the physical layer Deploying CFNC, a cooperative network 
with A'^s' sources offers throughput as high as 1/2 symbol per 
source per time slot (sym/S/TS), relative to 1/(2A^5) sym/S/TS 
in traditional cooperative system and l/{Ns + 1) sym/S/TS 
in GFNC-based cooperative system. Hence, CFNC has been 
applied to diverse cooperative networks [4]-[7]. 

2) Two-path relay. Relays usually work in the half-duplex 
mode, where a transmission process often requests two TSs for 
the relay to receive the signal and then forward it respectively, 
which results in a loss in spectral efficiency. To overcome it, 
in [8] and [9], a spectral efficient cooperative scheme named 
two-path successive relay was proposed. As illustrated in Fig. 
I, at each TS, the source S sends its symbol to one of 
the relays and the other relay forwards the symbol received 
from S at the previous TS. Therefore, only {L + 1) TSs 
are required to transmit L symbols from S to D. However, 
two-path relay scheme results in inter-relay interference (IRI), 
which is marked as the red lines in Fig. 1. As long as the IRI 
can be canceled, we can double the transmission rate without 
much sacrifice of the performance in terms of symbol error 
probability (SEP). In [10], network coding at relay nodes was 
applied to achieve inter-relay interference cancellation (IRIC). 
In [11], full interference cancellation (FIC) was proposed, 
which requires complicated detection. The orthogonality of 
the real part and the imaginary part of the symbol was utilized 
to distinguish the information and interference in [12]. More 
research on two-path relay in recent years can be found in 
[13][14]. 

Both of the methods above have improved bandwidth ef- 
ficiency and throughput of relay-based communication. Net- 
work coding deals with the increase of network size while two- 
path relay overcomes the half-duplex nature of relay nodes. 



Although CFNC allows simultaneous transmissions from all 
sources, 2L TSs are required to complete the transmission 
of L symbols for each source due to the half-duplex nature 
of relay. On the other hand, existing researches on two-path 
relay only consider single-source models. With the increase of 
network size, the single-source models in [10]-[14] cannot be 
directly generalized to multi-source systems with simultaneous 
transmissions at sources. 

Therefore, in order to further improve bandwidth efficiency, 
we propose a novel transmission scheme based on two-path 
relay in multi-source wireless networks. Our scheme combines 
network coding and two-path relay together in a (multi- 
source,2-relay,l -destination) system. To allow simultaneous 
transmissions of all sources, CFNC is operated before trans- 
mission. With the help of two-path relay, full-rate transmission 
of multi-source system is achieved. We use physical-layer 
network coding at relay nodes to deal with the IRI caused by 
two-path relay and the interference is successfully canceled at 
destination. Only L + 1 TSs are needed to send L symbols for 
each source. Compared with CFNC, our proposed cooperative 
transmission in (A^5,2,l) system achieves ideal throughput 
as high as 1 sym/S/TS, twice as high as that of CFNC in 
[3]. Compared with [10]-[14], multi-source model, rather than 
single-source, is introduced to two-path relay system without 
requirements for additional TSs and resource blocks. Besides, 
theoretical analysis and simulation results also demonstrate the 
improvement of our scheme's performance. 

The rest of this paper is organized as follows: Section II 
depicts the system model. Section III details the proposed 
transmission scheme and Section IV illustrates the IRIC algo- 
rithm. The theoretical analysis of the SEP and the throughput 
are presented in Section V. In Section VI, simulation results 
are given to verify the improvement of throughput in different 
scenarios. Finally, Section VII summarizes the paper. 

II. System Model 

As illustrated in Fig. 2, we consider a two-path relay coop- 
erative system that consists of multi-sources (5i, 5*2, ..., Sns)^ 
two relays (i?i,i?2) and one destination {D). One of the ap- 
plications is being adopted by the uplinks in cellular networks 
to promote transmission reliability for cell-edge users, shown 
in Fig. 2(b). Apart from extending coverage and achieving 
diversity gain, relay-based cooperative scheme could also 
reduce the required transmitted power at sources, which leads 
to lower inter-cell interference. Due to lower transmitted power 
at sources, we consider that the direct links between sources 
and destination (S-D) do not exist. All sources broadcast their 
own modulated symbols on the same RB in the same TS. 
Two relay nodes, working in the half-duplex mode, take turns 
to detect-and-forward (DF) received symbols to destination. 
Relay nodes also employ the same RB, which results in the 
IRI. Channel coefficients between 5, and Rm, Ri and i?2, 
R,n and D are denoted by hs^R,^, /112 and hji^^o where 
i = l,2,...,iV5 and m = 1,2. According to [15], when the 
coherence time of the channel is much greater than the symbol 
period, channel reciprocity (CR) can theoretically be assumed 
for channels in which uplink and downlink transmissions share 
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(a). Proposed tnulti-source cooperative system (b). Application in cellular network 

Fig. 2. System model 





TABLE I 




Notations 


Notations 


Description 


Ns 


the number of sources 


S^ 


the jth source, i = 1,2, ..., Ns 


Rm 


the mth relay, m = 1,2 


D 


the destination node 


^link 


channel coefficient, including: hs-R^ (S-R), hi2 




(R-B), hR^D (R-D) 


W 


AWGN, including: wr^ (received by R), wjy (re- 




ceived by D) 


(•)^ 


transpose 


e^ 


[01,02, ...,6ns]' precoding vector adopted by 




CFNC 


X 


[xi, X2, ...jXjVg]^, sources' original symbols 


Xr 


lxi^,X2^, ■■■,Xffg J*^, symbols sent by relay before 
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multiplied by 0^ 

the symbol sent by source % at time slot n 

the symbol sent by relay at time slot n coiTesponding 

to source i 

the estimation of x 

the symbol received by Rm and D respectively at 
TS n 

the PNC mapping function 
the real part and imaginary part of x 
the cardinality of modulation symbol set 
unit imaginary number 
0^x, the symbol from sources after CFNC 
X + X,., the superposition of desired symbol x and 
interference Xr 

0^z, the symbol received by R. in AWGN channels 
regardless of noise 

Hs^x+/ii2Xr, the superposition of desired symbol 
HgflX and interference h\2'^r 
0^z, the symbol received by R. in flat fading 
channels regardless of noise 
the modulation symbol set of x, s, 2, ■ ■ ■ 
the circular symmetiic complex Gaussian distribution 
with zero mean and variance cP' 
l/\/27r J°° exp(— i^/2)di, the Gaussian tail func- 
tion 

probability of error 
the mathematical expectation 
the noise power spectral density 
throughput 
the phase of h 



the same frequency spectrum. Therefore, regardless of relay 
nodes' mobility, we assume E[|/ii2p] = E[|/i2ip] = cr^2- 
The corresponding channel noise wh„, and wd are assumed 
to be additive complex white Gaussian noise (AWGN) with 
mean zero and variance A^o- All channel coefficients and 
noise are supposed to remain constant within one TS and 



vary independently from TS to TS. Throughout this paper, we 
mainly consider two kinds of channels: ideal AWGN channels 
with all channel coefficients being 1 and flat slow fading 
channels. Under normal circumstances, each receiver knows 
the channel state information (CSI) {Rm knows hs^R^ and 
hi2, D knows hft^^o)- 

In addition, we introduce some primary notations to describe 
our system, shown in Table I. 
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III. The Proposed Transmission Scheme 

Under the assumption that Ri and i?2 are assigned to 
relay messages for Ns sources, each source has continuous 
L symbols required to send to destination. To guarantee 
that different sources' symbols can be separated, CFNC is 
adopted here to establish a one-to-one correspondence between 
x=: [a;i,a;2, ....xjq^]^ and 9-^x where 6"^ = [^1,^2, ■■■^(^Ns]- 
Before transmission, the source symbol Xi from Si is multi- 
plied by Oi, i = 1,2, ...,Ns- The design of 8"^ is supposed 
to satisfy 8^x ^ 0^x' if x y^ x', which allows simultaneous 
transmissions from multiple sources. Among the different 
choices for Q^, we take it to be any row of the Vandermonde 
matrix; 

1 5i ... Jp- 

1 S2 ... 



S^^ 



-1 
1 



1 s 



Ns 



where {(5„} 

LA'S _ 



Ns 



„j-K{4:U-l)/i2Ns) 



NsxNs 



if Ns = 2*^ and 
{^u}u=i = eJ^(6«-i)/(3Ws) if Ns ^ 3 X 2^. This design 
is related to the linear complex field (LCF) encoder given in 
[16]. 

For simpUcity, we firstly introduce (2-source,2-relay,l- 
destination) system: 

A. (2,2,1) system 

Without losing generality, we assume that at time slot n, 
symbols from different sources are received by Ri and R2 
simultaneously forwards the symbol received at previous TS. 
While at time slot n + 1, sources' symbols are sent to R2 and 
_Ri forwards the previously received symbol to D, as shown 
in Fig. 3. 

Time Slot n: Si transmits 9iXi and simultaneously ^2 



J") 



» 



transmits 6*2 ajj where x\ denotes the modulated digital 
symbol sent by source i at TS n (i = 1, 2 and n := 1, 2, ..., L). 
At the same TS, i?2 sends the operated signal received at the 
previous TS n — 1 to D. Therefore, the received signal at _Ri 
and D are respectively given by 



vr; 



(n) 

Vd 



^hs^R^OixY' 



ISafli 



924"^ 



hi2{eix^r}+e2x';i>) + WR,, 



>h 



hR, 



D[ 



hx^r^ 



0A1) 



WD- 



(1) 



(2) 



We use subscript r to denote that symbols are sent by relay 
nodes. Hence Oixi'' +62X2 ' is the symbol sent by i?2 at TS 
n after PNC, which will be illustrated in the next section. All 

in') in) 

the symbols x\ , x^ transmitted by sources and relays are 



Fig. 3. (2,2,1) system model 



drawn from a finite alphabet Ax with cardinality \Ax\, which 
is determined by modulation mode. Maximum likelihood (ML) 
detection is performed at R and D as 



m 






{x. 



(») ^(») 



)]i? = arg min 



'e-4. 



h{hsiRiX\"' + hi2Xll') - 02{hs2RiX^2 



Ivr,- 

^12X2 J) 



(3) 



(ii' 



(») a(")- 



£> = arg min \\y}^' -hR^oi 



^i4:' 



924:') 



(4) 



where 6*1 and 62 are used to identify the symbols from 5*1 and 
5*2 respectively. If there is no interference (xi and X2 ), 
relay node is able to recover Xi and X2 [3]. Therefore, due 



(») 



(n) 

■ x) , rather than 



to the IRI, the relay just determines x. 
(x^" , Xj ). Then, we adopt PNC to deal with the interference. 
The core idea of PNC is encoding two sources' symbols at 
relay and decoding at destination [1]. In our scheme, the 

(n) (n) 

combinations of desired symbols and interference (x^ +xl ) 
are encoded through function / at relay and separated through 
function g at D. Let / : {xa + Xb\xa,Xb € Ax} — > Ax 



denote the many-to-one mapping function. We set x\ 



(n+l) 



'™(") 



^(")^ 



„("+!) 



/(4")- 



a; 



jyx-i ~T x-i j and Xo — j \'^2 ~'~*^2 
for i?i to forward within the next TS. 



Time Slot 71+ 1: 5*1 transmits 9ixl 

;ce 

("+i) 



5*2 transmits 62X2 , which are received by R2. At the same 



TS, i?i forwards Oix[''^'^^'' 



nX2^ 



to be the symbols 

and simultaneously 

by i?2. At the same 

to D. The received 



signal at D are given by 

IR^D\'^\X\^^ '+02X)2^ ')+WD, 

and the ML detector at D can be expressed as 



(n+l) ; 

Vd = '^' 



ii[T'\xt'''')D 



- arg mm 

' eA^. 
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I ("+i) 
IVd 



(5) 



(6) 



l-Ri D \ 



3rxt^'^ 



hx^t'^) 



With the purpose of detecting the sources symbols: 
Xi ,X2\ we introduce another many-to-one function g : 
{Ax , Ax ) — > Ax ■ According to the symbols detected at TS 
n and 71 + 1, given by (|4]) and ([6]), sources symbols can be 
recovered by 



(xi' 



(") An) 






(5(4:\4r''),.9(4:\4:+'')) 









= (,9(.Ti':,/(.Tr+4:')),ff(4:\/(4"' 



+ 4 ' 



)))• 



(7) 



An) 



The original symbol x\ is recovered through g{xl , /(a; 



-,(") 



.(") 



x(")^ 



(n+1) 



(«) 



)). Depending on the y^j and yjj stored at D, all the 
sources' symbols can be recovered through ([7]). The key issue 
is that whether the function / and g are able to successfully 
cancel the inter-relay interference (IRI), which will be detailed 
in the next section. 



B. (Ns,2,l) system 

The transmission scheme can be generalized to the multi- 
source (A^5,2,l) networks where Ns > 2. 

Time Slot n: According to the previous definition of x and 
0^, the I/O relationships are 






(8) 
(9) 

where ^sR = diag{hs^R^,hs^R^, ....hs^^R^)- x,"-* = 
[x\' , X2 , ■■■: Xn ]^ denotes the symbols sent by relay at TS 
n before multiplied by 9^. Through ML detection at relay, 
the symbol required to be forwarded by relay within the next 
TS is generated by x["+') = [f{x["^+x^;i^)J{4"^+4"J), 

•■•J Jy^Ns ~^ "^Wsr '^ ■ 

Time Slot n + 1: Ns sources continue broadcast symbols 
and simultaneously Ri forwards S^xJ." ' to D [17]. The 
received signal at D is given by 



(n+l) _ , oT^(n+l) I 

Vd ~ riRxD'o x;. + 



WD, 



(10) 



after which x}." can be obtained at D by: 



^(n+l) _ ^ ^^^ 



e^x 



\v'^'' 



h 0T„(n+l)||2 



Therefore, having x^ 



("+i) _ [^("+1) ^(«+i) 



.,x 



(n+l)lT 



together with xj- ' that estimated at last TS, each original 

[xi , X2 , ■■■, Xn ] ^ is recovered 



symbol x^"' in x 



r(") 



through 

Xi =g{xi',xi ') 



'a(") 



J") 



x(") 



gix\:>J{xr +£):>)) {^ = l,2,...,Ns). 



Finally, Ns sources successfully send L x Ns symbols 
during L + 1 TSs. As long as the destination stores current 
and former estimated symbols (xj." , xj." ), sources symbols 
x^") can be recovered through (fT2|) . Such transmission scheme 
achieves throughput L/{L + 1) sym/S/TS after the IRI cancel- 
lation (IRIC). When L is large enough, the ideal throughput 
reaches 1 sym/S/TS. 

IV. IRI Cancellation with PNC 

Based on the discussion in the previous section, the symbols 
received by relay from sources are interfered with the signal 
from the other relay. For the fact that relays are unable to ex- 
tract the desired symbols of sources, we introduce PNC to deal 
with the interference. We use function / at relay nodes and 
function g at destination to achieve IRI cancellation (IRIC). 
Since xl" = gix^^ f{Xi +Xj )), the most straightforward 
method is making use of the fact that 6 ffi a ffi 6 = a where © 



(n) 

denotes bitwise exclusive OR (XOR) operation. Hence x] is 
recovered through 

g{S^J(S^'>+x^y)=x^®fix^^+x^) 



®{xr®s:r) 



;,(") 



(13) 



To clearly illustrate the IRIC scheme, we make it simplify 
that hiink = 1- Without losing generality, all sources and 
relays adopt the same modulation mode: x^ , xf' £ 

Ax{n = 1,2,...,L). However, the symbols transmitted are 
0^x = s, rather than x. Satisfying the request for one-to-one 
mapping between x and s, Ns sources generate M^^ points 
in the constellation of s where M = \Ax\- Linear growth of 
Ns leads to exponential growth of cardinality of s. Therefore, 
we mainly consider modulation modes with low M such as 
BPSK(A/ = 2) and QPSK(A/ == 4), in order to achieve better 
anti-noise performance in multi-source system. 

A. BPSK modulation 

BPSK modulation generates two points in the constellation 
of each original symbol x. Regardless of high frequency 
carrier, the baseband digital symbol can be expressed as 
Ax — {1,-1} and M — \Ax\ = 2. Under the assumption 
that hiink = 1, the symbol detected by relay at TS n is 
Xj + x^' £ {2, 0, —2}. The estimated symbol is given by 



J") , ^(") 



[X^ ~rXi^ , X2 

= arg 






Tf^(n) I 



II v'^^ r\l / I 



r(n)) 



(14) 



„(") 



,(")|^(") ^(") 



where Az denotes {z = xl ' + xl^ '\xl \x\ 
we can express the function / as 



e Ax]- Hence 



/o 



x(") 



'Sn)\ 



x^r^®xf^ =\-\&^ ^if\ 



(15) 



from which /(2) = /(-2) = -1 and /(O) = 1. As long as 
satisfying a = fe © a © &, note that the so-called © not only 
means 1©1=:0, 1©0=1, but also varies based on different 
I/O relationships. 

At TS n + 1, 6^x1"+^) = e^[/(il"^ + il"^),/(x^"^ + 
(12) 4?'), -, !{x-nI + x~n\ )\^ is relayed to D. After detecting 
Xr , together with xj." stored at TS ri, another XOR 
operation is performed through p^ . In order to recover x^" , 
5(1, 1) = g(-l, -1) = -1 and g(l, -1) = g(l, -1) = 1 are 
supposed to be satisfied. Therefore, function g can be defined 
as 



x(«) 



x(«) 



^(") 



(") 



->) , .(n)^ 



gix^jixf^ + CO) = 1 - K' + m" + coi- (16) 

It is clear from (jTS]) that IRIC is completely successful 
and original symbols from sources are recovered at destina- 
tion. The difference between multi-source system and single- 
source system is the dimensions of vector x, among which 
pairwise dimensional components are relative independent 
and separated by Qj^ ^^ corresponding to the Ns sources. 
Therefore, without losing generality, taking 2-source system 
as an example, a complete transmission process is presented 
in Table II. The symbols in the last column are totally the same 
as the original symbols listed in the second column. Therefore, 
IRIC is completely achieved at destination. 



TABLE II 
IRICforBPSK 







TS 


n 


TS n + 1 






Vd 


Sources' symbols 


Received by R with IRI 


/(*!:^ +*!"') 


Vd 


ff(4"'.4"^'') 






x("' 


x(") 


xj,"' + xf") 


x;"+^' 


<v{"+l) 


x(") 






(1,-1) 


(1,1) 


(2,0) 


(-1,1) 


(-1,1) 


(1,1) 






(1,-1) 


(-1,-1) 


(0,-2) 


(1,-1) 


(1,-1) 


(-1,-1) 










TABLE III 
IRICforQPSK 








TS n 


TS n+l 


J") 
Vd 


Sources' symbols 


Received by R with IRI 


/(4"'+*.^"') 


yL"^^' 


5(4"^ 4"^'^) 


x'") 


x(") 


x("'+x("> 


x("+^) 




s(-) 


(l+j,l-j) 


(1+jM+j) 


(2+2j,2) 


(-i-i,-i+j) 


(-i-i,-i+i) 


(i-^i,i-Hj) 


(l+j,l-j) 


(-i+j,-i-i) 


(2i,-2i) 


(i-j,i-j) 


(i-i,i-i) 


(-i+j,-i-i) 


(i+i,i-i) 


(-i-i,i-i) 


(0,2-2j) 


(i+j,-i-j) 


(i-^j,-i-i) 


(-i-j,i-j) 


(-i-j,i-j) 


(-i-j,-i+j) 


(-2-2j,0) 


i-l-jA+j) 


(-i-i,i-^i) 


(-i-i,-i-^j) 


(-i-i,-i+j) 


{-\+jA+j) 


(-2,-2+2j) 


(-i+j,-H) 


(-i-^i,-i-i) 


(-i+j,-i+j) 



B. QPSK modulation 

QPSK modulation: Ax ^ {I + J, -1 + j, 1 - j, -1 - j] 
with Af = \-A.x\ = 4. QPSK symbols can be considered 
as the combination of two BPSK symbols: the in-phase and 
the quadrature-phase components. Estimated symbols 



^(") 



^(")^£(") 



(fft)) where yl^ 
2j 



x(") 






~(") 



)ii' are obtained through 



{2,-2,0,2j,-2j,2 + 2j,2-2j, -2 



2 — 2j}. Although Xj and x^" are not separable, the 



in-phase component 5R{: 
component 3{x 



x(") 



x(") 



(«) 



} and the quadrature-phase 
:" } can be easily separated. Therefore, 



it is not necessary to change the definition of / and g employed 
in BPSK modulation. The PNC operation at relay can be 
expressed as 

x^^'^ = fmt'^ +x^^})+jy. imt'^ +x^^}). (17) 

At TS n + 1, after detection of x^ , together with Xr 
stored at TS n, another XOR operation on divided in-phase 
and quadrature-phase components of the received symbols is 
performed through 



i;(")=g(5R{x(")},/Wi;l")+xi")})) 






(n) 



^(")l 



+ ]-Kg{'i{x'[^>}jmxT'+xTJ})). 



(18) 



Table III shows the transmission process of 2-source system 
with QPSK modulation. Similar to that of BPSK modulation, 
applying the same / and g at relay and destination respectively, 
the symbols from sources are completely recovered and entire 
transmission process is IRI free. 

C. General PNC process 

In practice, higher order modulation modes might be used, 
such as 16-QAM, 64-QAM [18] and so on. Apart from BPSK 
and QPSK modulation, our scheme can be extended to other 



modulation methods with different definition of / and g to 
achieve PNC mapping. To summarize, a general PNC process 
consists of the following steps: 

1) At each relay, for each source Si, many-to-one modulation 
mappmg, j: {x\ ' ^ x\ ') -^ x\ 



92x("+^) + 



2) Form the signal for relay to send: d\x]^ ' 

, a T,("+i) 

3) At D, for each source Si, one-to-one demodulation map- 

/ (ra+l) {n)-. ^ (n) 

ping, 9- (< \A^)^A ■ 

The designs of / and g vary based on different modulation 

modes, which are detailed in [1]. 

Now we consider a more practical system under flat slow 
fading channels. We assume that CSI is available at receiver 
nodes. At TS n, relay node i?i is aware of the CSI of 
link S-R:{hsiRi,hs2Ri,---,hst,^Ri) and the CSI of link 
R-R:hi2. Hence, Q can be generalized to A^5-sources sys- 
tem, expressed as 



m 






),(x. 



^(") ^(") 



(") ^(") 



Ns ' ^Ns 



arg 



mm 



e-4^ 



I (") 



2r h •■•' {^ 

- e^(H5Kx(") 



)]r = 
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We still mainly consider x\ + i,- . The symbol correspond- 
ing to source i for relay to forward can be obtained through 
fi^i + ^i )■ Ths IRIC at D is the same as the algorithm 
discussed above. In general, the total transmission process is 
similar to that in the system free from channel fading. 

V. Lower Bound for System Throughput 

In this section, we mainly analyze the throughput perfor- 
mance for the (A^5,2,l) system as described above. With 
the purpose of giving proof to verify that the proposed 
transmission scheme with IRIC can achieve 1 sym/S/TS under 
ideal circumstance, we assume that hunk = 1- All the additive 



complex White Gaussian noise (AWGN) w -- C7V(0,A'^o)- 
Both sources and relays adopt the same modulation mode, 
i.e., Ax with M = \Ax\- 

A. SEP of R-D link 

Firstly, we consider the SEP performance of R-D link. In 
traditional transmission procedure with single-source, the SEP 
of source symbol Xa{xa G Ax) can be calculated as 



Pe{Xa} = ^ P{Xa "> Xb} , 






(20) 



where P{xa — > Xf,} denotes the pairwise error probability 
of mistaking x^ for Xa. However, deploying CFNC in multi- 
source system, the transmitted symbol is 0-^x = s, which 
contains the information of all sources. Ns sources generate 
M^^ points in the constellation of s. However, the incorrect 
detection of s does not represent the error of all sources' 
symbols [I'l, a-2, ...jXjVg]. As an example to explain the idea, 
we assume Ns ^1 and BPSK [M ^ \Ax\ ^ 2) modulation 
is adopted by senders. Let 6i ~ 1 and 62 = j, as illustrated 
in Fig. 4. M^^ = 4 points exist in the constellation of s. 
Without losing generality, we assume the sources symbols 
are [—1,1], and thus s = 0i x (—1) + ^2 x 1, marked in 
the constellation. Even if the detected s ^ s, for example 
s = 6*1 X 1 + ^2 X 1, equivalent to x = [1, 1], from which one of 
the two sources' symbol is successfully recovered. Let random 
variable f denote the number of symbols failed to be recovered 
among Ns sources. Therefore, the probability distribution of 
C (0 ^ C ^ ^s) is given in Table IV where Pe denotes the 
SEP of Sa- The probability Pe-k [k = 1, 2, ..., Ns), denoting 
k of Ns symbols failed to be recovered at receiver, is given 
by 



P. 






Sb\ ||Xa -Xfa 



lo = fc}, 



(21) 



where Sa = &'^Xa,si, = 0-^Xf,. ||Xq — Xb||o indicates the 
zero-norm which denotes the number of nonzero elements of 
Xq — X(,. Hence, the SEP of R-D link for each source can be 
expressed as 

, Ns 
P.-RD - E{^}/Ns - TT I] fc X -Pe-fc, (22) 
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where E{^} is the mathematical expectation of ^. For any 
given constellation, the SEP in AWGN channel can be esti- 
mated by union bound[19]: 
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(23) 



where iVo is the noise power spectral density of the channel 
and (T^il denotes the Euclidean distance between s^ and Sb in 
the constellation. 

We assume the uniform distribution of s due to one-to- 
one mapping between s and x, i.e. P{sa} = l/l-^sl (a = 




X s = e' X 

Fig. 4. The constellations of symbols before and after CFNC 

TABLE IV 
Probability distribution of g 
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1,2,..., I A I). Based on ([23]), Pe-k is given by 
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according to which (|22|) can be rewritten as 
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B. SEP of S-R link 

Secondly, we consider the SEP of S-R link. Without losing 
generality, we set Nq-sr ~ ^o-i?,D = .^o- Links are free 
from channel fading and the signal received at R is the 
superposition of each node's symbol. The transmissions from 
Ns sources and one of the relays to the other relay can be 
treated as (1 -source, 1 -relay) link where the symbol sent by the 
equivalent source is 0^(x + Xr). Therefore, the analysis for 
R-D link can be applied to the derivation of SEP for S-R 
link. However, one of the differences between S-R and R-D 
link is that the symbol is 0^(x+Xr) = 6^z, rather than Q^x. 
Let y = Q^z. The most significant difference between s and 
y is that any Sa G -4^ appears with equal prior probability 
due to one-to-one mapping from x on the condition that 
each symbol in Ax has the same prior probability. However, 
the prior probability of ya G Ay varies due to many-to-one 
mapping (x, x^) — > (x + x^). For instance, deploying BPSK 
modulation, Az = {—2, 0, 2} with corresponding probability 
of 0.25, 0.5, 0.25, respectively. Therefore, when calculating 
Pe-k, P^ can be rewritten as 



(26) 



l-4yl / l-4„| / (y) N \ 

\||z„-z,||o=fc / 



Av) 



where P{ya} is the priori probabiUty of ya and crj^f^' denotes Appendix for the proof) according to (|25p and (pS)) : 
the Euclidean distance between ya and yi,. Hence the SEP of 
S-R link is given by 
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C. Lower bound for throughput with CFNC 

Traditional CFNC [3] requires 2L TSs to transmit L con- 
tinuous symbols: at odd TSs, sources send symbols to relays; 
at even TSs, relay nodes forward the symbol to D while 
sources keep silence. Since we have Pe-SR and Pc-rd, the 
throughput of S-R-D with traditional CFNC can be estimated 
by 



1, 



TcFNC ^ ^[{l ~ Pe- 



SR 



){1 ~ Pe-RD) + P'] 



(sym/S/TS), 
(28) 

where throughput T is defined as the number of successfully 
transmitted symbols per source per TS (sym/S/TS) within one 
transmission cycle. When calculating P^-rd, Nq in (p5|) is 
supposed to be the value corresponding to the SNR after maxi- 
mal ratio combining (MRC) of SNRajD and SNRa^d. Due to 
no IRI in traditional cooperative communication with CFNC, 
-^e-SR ^'^^ ^^ obtained through (P5|) , following the analysis of 
link R-D, rather than ([27]) . The factor ^ means that two TSs 
are required to send one symbol from source to destination. 
The first term {l — Pl,_gj^){l — Pe-RD) represents no incorrect 
decision at both R and D, i.e. Xa — > Xa — J- Xa, while the 
second term P' denotes the probability that one symbol is suc- 
cessfully recovered at D while incorrect decisions occur within 

.u . ■ • • .1- SR RD , , , 

the transmission, i.e. the case Xq ^- a-;, — > Xa [Xa ^ xi,)- 
Therefore, P' is given by 
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D. Lower bound for throughput with proposed scheme 

Then, we consider the throughput of the proposed trans- 
mission scheme with IRIC. According to (H^), the detection 
of original symbols depends on current signal received at D 
and the signal arrived at former TS as well. Besides, at the 
first TS, there is no IRI at i?i since R2 keeps silence, which is 
equivalent to traditional CFNC scheme. Hence, the SEP of the 
symbols transmitted at TS 71 (1 < n < L) can be expressed 
as 

pin) \ 1 - (1 - ^e-5fl,)(l - Pe-RD) - P' n=\ 

"^" 1 1 - (1 - Pe-SR){^ - Pe-RD? - P" ne [2, L] 

(31) 

i + 1 TSs are required to transmit L symbols for each source. 
Therefore, the throughput of proposed transmission scheme 
with IRIC is given by 
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where the first term (1 — Pe_5/j)(l — Pe-Ru)^ indicates that 
the respective symbol is successfully transmitted from sender 



S^-^R 



TS n 

to receiver without errors in each link: R — > D, 
and R — > D. P" denotes the probability that incorrect 
decisions occur at least twice among all three links while the 
source's symbol still can be recovered by D. The calculation 
of P" is divided into two cases: 

7 1 7 1VT ■ r> TS n r^ • SR RD 

1) case 1: No error occurs m R — >^ D, i.e. Xq — > a:;, — > 
Xa {xa ^ Xb) at TS n + 1. Similar to the analysis on traditional 
CFNC, the probability of case 1 can be expressed as 
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where A.w = {st\sb e As, {yia)N ^ Mn} and {Xa)N P'" ^TT 51 E ^^^^^^ E ^i^" ^ yb]P{sb ^ Sa} 
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denotes the A^th element of x^ 

Under the assumption that A^o- 
have Pe-RD 

based scheme, which results in Tcfnc ~ \[{^~ Pe-RoY' + 
P']. Hence, we achieve the lower bound for throughput (see 
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where Aa^N == {yblvb e Ay, {za)N ^ (zb)^}- 
2) case 2: Transmission errors occur in R 
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D. We 
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assume that x\r , the A^th element of x^ 
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is sentenced to x^ in R-D Unk at TS n. In order to recover 

(n) 

x\r simultaneously sent by sources at TS n, the symbol 
x^" detected by Z? at TS n + 1 should satisfy 

where x]^ and x^ are the correct symbols supposed to be 
detected by D without errors. Hence, the probability of case 
2 can be expressed as 
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IT 

yb — Sa + Sb and s'^. is the correspondence to 

4, = [■••:4„,.-]'^' which satisfies xn = 5««,.,4„J- 
Replacing P{- • • } with Q{- ■ ■), Pcase-2 is given by 
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According to (lp]),(1271),(I25]) and ^, we define Pe^Ro, 
Pe-sR, P'" and Pcase-2 as the bounds for Pc-rd, Pe-SR, 
P'" and Pcase-2 respectively by replacing P{---} with 
g(- • • ). Besides, F™«e-i = {1 - Pc-rd) ■ P"' ■ Let 



-PrrvK/^-l — (1 — Pe-BD) ■ -P • 



case— 1 

case—1 



(39) 



Based on the fact that P" = Pcase-i + Pcase-2, we obtain 
the lower bound for r„e«,: 

Tnew > (1 — ^e-Sfl,)(l ^ Pe-Ro) + ^case-1 + ^case-2, 

(40) 

and the proof of ">" is similar to the derivation of CFNC- 
based scheme. 

In AWGN channels without fading, the union bound for 
SEP, obtained from (P^ . approaches the actual value. There- 
fore, the lower bounds can be the approximations of practical 
throughput, especially in high SNR regime. When SNR tends 
to infinity (i.e. A^o tends to zero), the SEPs tend to zero and 
we can achieve 

TcFNC = Y^va\[{l-P'c_sR){l-Pe-RD) + P'] = \ (41) 
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from which the throughput of 1 sym/S/TS is verified. 



E. Practical system with Rayleigh fading channels 

Now we consider a more practical system where the chan- 
nels are assumed to be flat Rayleigh fading. We have hs^R^ ^ 
CAf{0,afJ,hi2,h2i ^ CAf{0,al^),hR^D ^ CAA(0,af^), 
{i ^ 1,2,..., Ns,m^ 1,2). 

We still firstly analyze the link R-D, where the transmitted 
symbol is Hr^dQ^x,-. Each vector's length in the constel- 
lation of s is multiplied by IHr^dI which follows Rayleigh 
distribution. Therefore, P{sa — > si,\hR^]:)} can be expressed 
as 



P{Sa -^ SblhR^o} < Q 
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where d^^; denotes the Euclidean distance between Sa and st 
in AWGN channel. According to Chernoff bound [20], we 
acquire 
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Probability density function (pdf) of Rayleigh distribution is 
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Besides, the phase of hR^jj follows uniform distribution on 
intervals [0, 27r] and is independent of the distribution of 
I^r^dI Hence, let v,nD = l^fl^ul and ZHr^d = a, we 
can acquire 
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with 



2JVo 



in (1251), the SEP 



Replacing Q , ^.-^ 

bound for Rm-D link is attainable. Pr-rd ~ ^{Pe-RiD + 
Pe-R-iD) owing to the equivalence of two relays. 

Then, we consider the performance of S-R link. Regardless 
of noise, the signal received by relay is 0-^H5j7X + 0^/ii2Xr. 
Let z = HsflX + /ii2X,. and y = Q^i. Then, (P7)) can be 
rewritten as 
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where djf^ denotes the Euclidean distance between ija and yi,. 



Since all channels are independent of each other, we have 
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where H^i? = diag{hsiB.„,,hs2R^, 

|/iS,fl„| = 'yim,|/ll2| = Wl2 {i = 1,2,. 

The phases of channel coefficients are Zhs^R. 
Zhi2 = ai2- In addition, 
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and 
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n and 






dT = |e^(Hsflx„ 



hi2Xa^ - UsR^b - /li2XhJ|^, (49) 



where HsRXa + /ii2Xa^ = Zo and HsRXb + ^i2Xfc^ = Zf 
Since all the Ui,„ and a^m are independent of each other, pS)) 
can be obtained through one by one integration. To reduce the 
complexity of integral, we make some simplifications, e.g., if 
the sources stay close during a period of time, we may claim 
the same channel coefficient for each of them. Besides, due to 
the large scale of constellation, we can only use d™*", rather 
than d^^ , for each ija where d™*" is the minimum Euclidean 
distance between ija and other points in the constellation. Since 
d™™ mainly determines the SEP performance of Pe{ya}, the 
mean error can be low enough. For any given Xq, Xa , xj, and 



x;,^, by replacing Q 



/2NZ 



with dm), the SEP bound for 



S-R link (|T7)) can be achieved by calculation. 

The subsequent analysis of throughput is the same as that 
in AWGN channels, based on Pe-SR and Pc-rd- The perfor- 
mance of our proposed scheme in Rayleigh fading channels 
through simulations is presented in the next section. 
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Fig. 5. Throughput for (2-2-1) system in AWGN channels 
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Fig. 6. Throughput for (3-2-1) system in AWGN channels 



VI. Simulation 

A. Simulation Setup 

We apply the transmission scheme to the uplinks in cel- 
lular network and the simulation parameters are configured 
according to the 3GPP LTE specifications [21] [22] and specific 
parameters related to such transmission scheme as well, which 
are shown in Table V. 

TABLE V 
Simulation parameters 



Parameters 



Value 



Modulation Mode 

Packet Size 
User Antenna Gain 

Noise Power 
Path Loss 



BPSK.QPSK 

128 bits 

OdBi 

LCF[14] 

-174 dBm/Hz 

15.3 + 37.61ogiod 



Uplink Receiver Type MRC(Maximum Ratio Combining) 

Two scenarios of channel are considered: AWGN channel 
and Rayleigh fading channel. The design of O^ is the same as 
the description in section III. Through ML detection, symbols 



are detected at both relay nodes and destination. Assuming that 
the noise power spectral density iVo remains unchanged during 

the entire transmission process, Q f" ] is equivalent to 

Q ( ^\/SNr) where SNR is the received SNR and d'^^ is 

the Euclidean distance between Sa and Sb in the constellation 
where E{||s|p = 1} (i.e. transmitted power is normalized to 
be 1). Given modulation mode, the number of the sources (Ns) 
and the design of 0^, the throughput is only determined by 
SNR. We still assume iVo-Sfl = Nq-rd = Nq. Therefore, 
we mainly consider the performance of SEP and throughput 
in various scenarios for practical SNR values. Additionally, 
since the linear growth of Ns leads to exponential growth of 
cardinality of Q^x and 6^(x + Xr), which results in deterio- 
ration of SEP and throughput, we mainly consider the systems 
with 2-source and 3-source. To verify the improvement of 
the performance and illustrate the difference between 2-source 
and 3-source, the throughput in this section is defined as the 
number of successfully transmitted symbols from all sources 
per TS (sym/TS), rather than sym/S/TS. 



10 



-2-2-1 (Proposed Scheme) 
-2-2-1 (CFNC) 
-3-2-1 (Proposed Scheme) 
-3-2-1 (CFNC) 




5 10 

SNR [dB] 



10" 




-2-2-1 (Precoder design with CFNC) 
-3-2-1 (Precoder design with CFNC) 
-2-2-1 (Proposed Scheme) 
-3-2-1 (Proposed Scheme) 



10 15 20 

SNR [dB] 



Fig. 7. Throughput perfonnance in Rayleigh fading channels 



Fig. 9. SEP comparison in AWGN channels 
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Fig. 8. Throughput comparison for various channel SNR settings 



Fig. 10. SEP comparison in Rayleigh fading channels 



B. Results and Discussions 

Firstly, we compare the performance between proposed 
scheme and traditional CFNC scheme. CFNC [3] needs two 
TSs to complete one transmission. In our proposed scheme, 
with two relays alternatively forwarding the symbols to D, 
only (L + 1) TSs are required to transmit L symbols. BPSK 
modulation is adopted for all nodes. We evaluate the through- 
put performance for (iV5',2,l) systems in AWGN channels. 
In Fig. 5 and Fig. 6, the average throughput of proposed 
scheme based on two-path relay with IRIC in (2,2,1) and 
(3,2,1) systems is compared with that of traditional CFNC. 
Even when SNR is as low as -5 dB, the throughput of new 
scheme is higher than the upper limit of traditional CFNC. 
With the increase of SNR, the throughput of proposed scheme 
approaches the upper limit, which is 2 sym/TS in 2-source 
system while 3 sym/TS in 3-source network. Furthermore, 
we calculate the theoretical values of throughput based on 
(pO)l and (|40)l . marked as Theory, and compare them with the 
simulation values under different SNR. To clearly investigate 
the tendency of the lower bound derived above, SNR varies 
from 5 dB to 30 dB due to large scale of constellation for 
3-source system. It can be observed that the theoretical values 



and the simulation values are almost identical in high SNR 
regime, which verifies the derivation in Section V. Moreover, 
the theoretical values demonstrate that even the lower bounds 
achieve Ns sym/TS, not to mention the actual values. Then, 
we consider a more practical system where the channels are 
assumed to be flat Rayleigh fading. Although the throughput 
decreases compared with that in AWGN channels due to 
channel fading, the new scheme can also reach the upper limit 
in high average SNR regions, as shown in Fig. 7. 

Secondly, throughput comparisons for various SNR settings 
in Rayleigh fading channels are given in Fig. 8. The relay 
nodes always locate between sources and destination, so link 
S-R or R-D achieves better link quality than that of direct 
link S-D. Link quality can be represented by SNR. Since 
R is located either close to sources, close to D or not 
close to both, we set corresponding SNRs (7sfl,7i?,D) in 
logarithmic-scale are (SNR+10 dB,SNR), (SNR,SNR-f 10 dB) 
and (SNR,SNR), respectively. According to the simulation 
results, the relays close to sources (SNR+10 dB,SNR) are 
supposed to have high-priority in relay selection, especially 
in high SNR regime. Furthermore, better link quality of S-R 
requests lower transmitted power at sources, which reduces 
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TS on the same RBs. In order to deal with the IRI caused by 
relays, PNC is applied at relay nodes and IRI is successfully 
canceled at destination. Theoretical analysis is presented to 
verify the performance of proposed scheme. Lower bounds 
for throughput in AWGN channels and Rayleigh fading chan- 
nels are given to estimate throughput of such system. The 
simulation results in various scenarios illustrate that this new 
scheme achieves higher throughput than that of other schemes 
in either AWGN or Rayleigh fading channels. 

Appendix 

The proof for ^■. 

Through ^, we define P' to be 



Fig. 11. Thi'oughput performance adopting QPSK modulation 



inter-cell interference caused by cell-edge users. 

A novel transmission scheme based on precoder design is 
proposed in [7], which also achieves symbol rate as high 
as 1 sym/S/TS regardless of noise and interference. During 
traditional transmission process, the symbol is transmitted one 
by one separately. In [7], source's symbol transmitted in chan- 
nel is the superposition of two original modulated symbols, 
which doubles the ideal throughput. Therefore, we evaluate 
the performance of symbol error probability (SEP) between 
our new scheme and the scheme proposed in [7], shown in 
Fig. 9 and Fig. 10. Although the scheme based on precoder 
design is free from IRI, half of the symbols are detected based 
on S-D link in [7]. However, the link quality of S-D is 
worse than that of S-R-D, which apparently deteriorates the 
whole performance. Therefore, in either AWGN or Rayleigh 
fading channels, it can be observed that our proposed scheme 
achieves better SEP performance than what is proposed in [7]. 

All the simulations above are based on BPSK modulation. 
To check modulation other than BPSK, we test the throughput 
performance of proposed scheme adopting QPSK modulation. 
From Fig. 11, although many transmission errors occur under 
low SNR due to large scale of the constellation, the throughput 
can also achieve the upper limit when SNR is large enough in 
either AWGN or Rayleigh fading channels. 

In general, the simulation results have verified that through- 
put of the proposed scheme in this paper is able to reach 
the upper limit (1 sym/S/TS). Compared with other method 
which can also attain 1 sym/S/TS in iVg— source system, our 
proposed scheme achieves better performance. 

VII. Conclusion 

In this paper, we propose a novel cooperative communica- 
tion scheme based on two-path relay with IRIC to improve 
throughput in multi-source wireless networks. The two-path 
successive relay scheme ensures that the source continuously 
sends symbol to the two relay nodes alternatively at any 
TS. Deploying CFNC, the single-source system with two-path 
relay is generalized to multi-source system, from which all the 
sources are able to broadcast their own symbols in the same 



P' 



Ns\As 



Ns 1-4.1 

EE E 

N=l a=l SbeAa 



As) 



Q 



(50) 



Therefore, we have 1 — Pc-rd > 1 ^ Pe-R.D and P' < P'. 
Due to the fact that a > b and c < d cannot directly lead 
to a + c > b + d, our purpose is to verify that Tcfnc ^ 

W-Pe-B.D?+P']- 

The throughput Tcfnc can be rewritten as 



Tcfnc — x[1 ^ 2Pe-B.D 



P: 



e-RD 



P'] 



= -[1 - Pe-RDil - Pe-RD) + [P' - Pe-Ro)]. 



where 



-Pe-RD{1 - Pc-Rd) > -Pc-Rd{1 - Pe-RD) 



(51) 



(52) 



Besides, let Pe-RD + A;,^ = Pe-RD. Pe-RD + ^'rd = 
P2_^^ and P' + Ap, = P' (0 < Ard, A'^^, Ap. < 1), we 
have 



P 



e-RD 



'J-Pe- 



RD 



e-RD 



A 



RD 



a-p 



e-RD) 



(53) 



RD- 



From ((52l) and ((53l) . we have Ard — A^^^ > 0. Moreover, 
note that Pe-RD consists of all the cases when transmission 
errors occur, i.e. Xa —>■ xi, —> x^ (xa 7^ xt,,xi, 7^ Xk), which 
contains P' : Xa -^ x^ — > Xa {xa 7^ xi,). Therefore, it is 
obvious that A^^^ > Ap/, so is Aro — Api > 0. Hence, the 
second term of (ISTI): 



P' - P, 



e-RD 



= P'- Pe 



-e-RD + Ard - Ap, >P'- Pe-RD, 

(54) 

together with ([5^ . the lower bound for Tcfnc is obtained 
by 



Tcfnc — -^[(1 ^ Pe-RD) 



P'] >-[{!- Pe-RD)^+P'] ■ 
(55) 



The lower bound for Tcfnc in (PH)) is verified. Furthermore, 
the lower bound for Tnew can also be verified through similar 
method presented above. 
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